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Abstract 

We investigate the behaviour of the nuclear forces as a function of the light quark masses (or, 
equivalently, pion mass) in the framework of chiral effective field theory at next-to-leading order. 
The nucleon-nucleon force is described in terms of one and two-pion exchange and local short 
distance operators, which depend explicitly and implicitly on the quark masses. The pion propa- 
gator becomes Coulomb-like in the chiral limit and thus one has significant scattering in all partial 
waves. The pion-nucleon coupling depends implicitly on the quark masses and we find that it 
becomes stronger in the chiral limit. There is a further quark mass dependence in the S-wave four- 
nucleon couplings, which can be estimated by means of dimensional analysis. We find that nuclear 
physics in the chiral limit becomes natural. There are no new bound states, the deuteron binding 
energy is £?£ L = 9.6 ± 1.9 MeV, and the S-wave scattering lengths take values of a few fm, 

acL^^o) = -4.1 ± 1-6+0 4 fm and acL( 3 5*i) = 1-5 ± 0.4 3 fm. We also discuss the extrapolation 
to larger pion masses pertinent for the extraction of these quantities from lattice simulations. 
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1 Introduction 



The application of (chiral) effective field theory (EFT) to few-nucleon systems has become a subject 
of vigorous research interest in the last years. One successful method is based upon the original idea of 
Weinberg [1] to apply the standard methods of Chiral Perturbation Theory (CHPT) to the kernel of the 
corresponding integral equation (Lippmann-Schwinger equation in case of two nucleons and Faddeev- 
Yakubovsky equations for three and more nucleons), and consequently solving it numerically. The 
necessity of summing the interactions to an infinite order is caused by the nonperturbative aspects 
of the problem (shallow nuclear bound states, unnaturally large S-wave scattering lengths), quite in 
contrast to the purely perturbative treatment possible in the pion and pion-nucleon sectors. The 
complete analysis of two- and three-nucleon systems and, in addition, of the a-particle has been 
recently performed up to next-to-next -to-leading order (NNLO) in the chiral expansion [2, 3, 4, 5]. 
Going to NNLO turns out to be sufficient to describe accurately most of the low-energy observables. 
CHPT is an effective field theory based upon the approximate and spontaneously broken chiral sym- 
metry of QCD and allows for a model-independent and systematic way of calculating low-energy 
hadronic observables. This is achieved by a simultaneous expansion in small external momenta and 
quark (pion) masses around the chiral limit, in which the light quarks (and, as a consequence, also 
Goldstone bosons) are massless. Since this expansion is well defined for sufficiently small quark masses 
(m q ), hadronic properties at low energies are expected not to change strongly in the limit m q — > 0. 
This feature is of crucial importance for the applicability of CHPT and is certainly valid in the purely 
pionic sector. In the chiral limit the interaction between pions becomes arbitrarily weak for vanishing 
external momenta. The complete one-loop analysis for the two flavor case has been performed by 
Gasser and Leutwyler in 1984 [6]. By now, the two-loop analysis of the irir scattering and other selected 
S-matrix elements and transition currents has become available, see e.g. [7]. In the pion-nucleon case 
the situation is similar, since the coupling of the Goldstone bosons to matter fields vanishes at small 
momentum transfer in the chiral limit as well. The loop expansion in the presence of the additional 
mass scale (i.e. the nucleon mass in the chiral limit) is somewhat more tricky, but various methods 
have been invented to allow for a consistent and symmetry preserving power counting. The purpose of 
the present work is to consider the nucleon-nucleon (NN) interaction in the chiral limit. Clearly, this 
case requires additional care because of the nonperturbative aspects of the problem and also due to the 
fact that the interaction between nucleons does not become weak in the chiral limit.* 4 In principle, 
at very low energies one may consider a theory of nucleons only interacting via contact interactions, 
the so-called pionless EFT (see e.g. [8] and references therein). While such a scheme has been shown 
to be a precise calculational tool, the link to the chiral symmetry is lost and we thus eschew such 
type of EFT here.* 5 We stress that the question about the M^-dependence of nuclear forces is not 
only of academic interest, but also of practical use since such studies of quark mass dependences can 
be used to interpolate lattice gauge theory data, which are usually obtained for quark (pion) masses 
much larger than physical ones. Furthermore, the deuteron and other light nuclei can be considered as 
laboratories to test the recent ideas of the time variation of certain fundamental couplings constants 
(see e.g. [9] and references therein). 

The first step in extrapolating the properties of the NN system for vanishing quark mass (exact chiral 
symmetry of QCD) was done by Bulgac et al. [10] a few years ago. They stressed that the vanishing 
pion mass may potentially lead to strong effects in the nucleon-nucleon P- and higher partial waves. 
One could even not exclude a priori the presence of bound states in some of the P-waves in the chiral 
limit. This indicates that in case of few-nucleon systems, the expansion around the chiral limit might 
be complicated or even not well defined. Based upon the one-pion exchange potential (OPEP) and 

* 4 Moreover, it even gets stronger in most channels, as we will show below. 

#5 In the pionless EFT the pion mass is considered to be a large scale and the limit M n — » cannot be performed. 
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taking into account only its explicit dependence on the pion mass, the authors of ref. [10] came to the 
conclusion that no bound states appear in P- and higher partial waves in the chiral limit. Performing 
first-order perturbation theory and using the deuteron wave-functions from some phenomenological 
potential models, Bulgac et al. found a reduction of about 1 MeV (that is by 50%) in the deuteron 
binding energy. The main conclusion of ref. [10] was therefore that ". . . physics of nuclei in a universe 
where M n = would be similar to what is actually observed." However, the interpretation of the 
results found in [10] requires some caution not only due to the fact that first-order perturbative 
results can not be trusted in that case,* 6 but also because only the explicit dependence on the pion 
mass has been taken into account. In the language of chiral effective field theory, taking into account 
only the explicit M^-dependence in the OPEP corresponds (roughly) to the LO approximation, which 
is not accurate enough to lead to any quantitative conclusions. 

More recently, Beane et al. re-analyzed the situation in the 1 S'o and 3 Si — 3 Z?i channels including 
also some implicit dependence of the NN interaction on the pion mass [11]. They considered the 
OPEP as an approximation to the NN interaction for distances larger than some matching radius 
For shorter distances the potential was approximated by a square well, which can be viewed as 
a smeared out delta-function counter term. In the 3 Si~ 3 Di-channel they also included an additional 
energy-dependent contribution in the short-range part of the potential. Adjusting the two parameters 
related to the short range part of the interaction, Beane et al. performed exact renormalization of 
the corresponding Schrodinger equation, i.e. the calculated low-energy observables do not depend on 
the matching radius R±. This way of handling the diverging Schrodinger equation is closely related to 
their previous work [12]. The situation in the chiral limit has been analyzed in two ways: (i) Keeping 
only the explicit M n dependence of the OPEP an increase of the deuteron binding energy by about a 
factor of two has been found, ~ 4.2 MeV. (ii) The authors of ref. [11] also made an attempt to take 
into account the implicit M^-dependence, which changes the strength of the OPEP in the chiral limit. 
They further assumed that these corrections can be represented by the leading logarithmic terms in 
the chiral expansion.* 7 Quite surprisingly, the results for deuteron binding energy were found to 
change completely after taking into account the implicit M^-dependence in this manner. It turned 
out that deuteron becomes unbound for the pion mass smaller than ~ 100 MeV. Another interesting 
issue discussed in [11] refers to the possibility of matching the results obtained in chiral effective field 
theory with the ones from lattice QCD [13], which requires extrapolation away from the physical value 
ofM*. 

While this paper was in preparation, an NLO analysis of the quark mass dependence in the 1 Sq and 
3 Si~ 3 Di channels by Beane and Savage based upon the power counting introduced in ref. [11] has 
become available [14], improving and extending the results of ref. [11]. Apart from the OPE and 
short-range terms in the effective potential, the leading TPE contribution (in the chiral limit) has 
been taken into account. The authors found evidence for nearly all possible scenarios such as a bound 
or unbound deuteron and di-neutron in the 1 5o channel and could not make any predictions for the 
corresponding scattering lengths due to the lack of knowledge of the low-energy constants related 
to contact terms with one insertion proportional to M%. We will come back to this work and to 
comparison to our results in the next-to-last section. 

The aim of this work is to perform a complete NLO calculation in the framework of a modified 
Weinberg power counting and to try to clarify the presently somewhat controversial situation about 
the chiral limit of the nucleon-nucleon interaction # 8 . To do that, we have to consider apart from the 

# 6 As pointed out in [11], one gets a deuteron binding energy Bd ~ 4.1 MeV if the pion mass is set to zero in the AV18 
potential. Thus the effect is of the opposite sign to the one found in [10] using first-order perturbation theory. 

* 7 Note that this assumption goes beyond standard CHPT and is not applicable in a general case. We will come back 
to this point later on. 

#8 This is similar to what was done in [14] but differs markedly in some aspects from that work as will be shown below. 
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OPEP and contact interactions also the leading two-pion exchange potential (TPEP). We will discuss 
in detail the renormalization of the NN potential at NLO since this is of utmost importance for such 
type of analysis. We will not use of the assumption of the dominance of chiral logarithms but rather 
work with the complete CHPT expressions at this order. The corresponding low-energy constants 
(LECs) will be taken from the analysis of various processes. We will demonstrate that using only 
the leading logarithms leads to results for the renormalization of the axial coupling qa incompatible 
with the present analyses of the ttN system and of the process irN — ► irnN. We also perform an 
extrapolation of the 1 Sq and 3 Si scattering lengths and deuteron binding energy Bjy for the values of 
the pion mass in the range < M n < 300 MeV. 

Our manuscript is organized as follows. In section 2 we consider the chiral effective nucleon-nucleon 
potential at next-to-leading order, with particular emphasis on the aspects of renormalization per- 
tinent to derive the desired quark mass dependences. Such an analysis is not yet available in the 
literature. We also demonstrate the equivalence between the S-matrix approach used to analyze the 
pion-nucleon sector and the non-covariant projection formalism employed to construct the NN po- 
tential. We systematically work out all explicit and implicit quark mass dependences of the various 
contributions due to one- and two-pion exchanges as well as contact interactions. Results for the NN 
phase shifts, the deuteron binding energy and the S-wave scattering lengths are collected in section 3. 
Will also discuss in detail the differences to the earlier EFT work [11, 14]. A brief summary and 
further discussion is presented in section 4. 

2 Nucleon-nucleon potential at next— to— leading order 

In this section, we remind the reader on the structure of the NN potential at NLO and perform its 
complete renormalization, which is of crucial importance for considering the chiral limit. This has 
not been done in our previous work [3], since we were only interested in constructing the potential 
for the physically relevant case. Since different methods have been applied in the literature to define 
the effective interaction it appears necessary to begin our considerations with a brief overview of some 
commonly used formalisms. This topic is also of technical importance because the pion-exchange is 
directly linked to the pion-nucleon sector, and thus consistency with the field theoretical NN potential 
has to be demonstrated. The reader more interested in the explicit form of the potential at NLO 
exhibiting explicit and implicit pion mass dependences may directly proceed to section 2.4. 

2.1 Construction of the NN potential from field theory 

The construction of a potential from field theory is a well known and intensively studied problem in 
nuclear physics. Historically, the important conceptual achievements in that direction have been done 
in the fifties in the context of the so called meson field theory. The problem can be formulated in 
the following way: given some field theoretical Lagrangian for interacting mesons and nucleons, how 
can one reduce the (infinite dimensional) equation of motion for mesons and nucleons to an effective 
Schrodinger equation for nucleonic degrees of freedom, which can then be solved by standard methods? 
It goes beyond the scope of this paper to discuss the whole variety of the different techniques which 
have been developed to construct an effective interactions, see ref. [15] for a comprehensive review. 
We will now briefly introduce two methods which will be important for our further considerations. 
The first method is closely related to the field theoretical S-matrix given (in the interaction represen- 
tation) by 

S a p = S a p-2^iS{E a -Ep)T aP , (2.1) 
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where the T-matrix T Q/ g satisfies: 

T a p = V aP + V ai {E p - £ 7 + ie)-% p . (2.2) 

Inverting the last equation one obtains the following equation for the effective potential V a p: 

V a/ 3 = T afS - T ai {E p - £ 7 + ie)-%fs + . . . , (2.3) 

where the ellipsis refer to higher order iterations in the T-matrix. The effective potential V a p can 
now be obtained in terms of a perturbative expansion for the T-matrix. For example, in the Yukawa 
theory one can calculate T a/ 3 (and also V a p) in terms of a power series expansion in the squared meson- 
nucleon coupling constant g 2 . It has to be pointed out that the potential V a p in eq. (2.3) is not defined 
unambiguously since only on the energy shell T-matrix elements are known and well-defined in field 
theory. In fact, at any fixed order of the perturbative expansion (2.3) one can perform an arbitrary 
off-the-energy shell extension of the potential, which will then affect even on the energy shell matrix 
elements at higher orders. In addition, one has a freedom to perform field redefinitions, which changes 
off shell S-matrix elements and, as a consequence, off shell matrix elements of the potential. Such an 
ambiguity of the effective potential is a general property and does not pose any problem, since the 
latter is not observable. Only observable quantities, i.e. on shell S-matrix elements, binding energies, 
. . ., are defined unambiguously in quantum field theory. 

The method introduced above has been used by Kaiser et al. to derive the NN potential from the 
chiral Lagrangian [16, 17] and will be refered to in what follows as the S-matrix method. Its close 
relation to the scattering amplitude makes it possible to apply the standard field theoretical techniques 
for calculation the effective interactions. 

Another well known scheme is often refered to as the Tamm-Dancoff method [18, 19]. The starting 
point is the time-independent Schrodinger equation 

(H + H I )\V)=E\V), (2.4) 

where |^) denotes an eigenstate of the Hamiltonian H with the eigenvalue E. One can now divide 
the full Fock space into the nucleonic subspace \<p) and the complementary one \ip) and rewrite the 
Schrodinger equation (2.4) as 

V Hr, r}H\\(\4>)\_ ( \<f>) \ 

XHr, XHX ) \ \ip) ) ~ \ m J ' [ b) 

where we introduced the projection operators r] and A such that \<f>) = r/|^), \tp) = \\^)- Expressing 
the state \tp) from the second line of the matrix equation (2.5) as 

i* = m^m™ ■ < 2 ' 6 > 

and substituting this into the first line we obtain the Schrodinger-like equation for the projected state 

(H + V cS (E))\<i ) )=E\<i ) ), (2.7) 
with an effective potential V e g (E) given by 

V cS (E) = nHjr) + VH I X E _ 1 XHX \H IV . (2.8) 

Eq. (2.7) differs from the ordinary Schrodinger equation by the fact that the effective potential V c g (E) 
depends explicitly on the energy, as it is obvious from the definition (2.8). In fact, such an explicit 
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energy dependence can be eliminated from the potential as discussed e.g. in ref.[21]. Indeed, the 
effective potential V e s (E) up to the fourth order in the coupling constant g in the Yukawa theory with 
Hj = gH l is given in a symbolic form as 

(V eS (E)) mn - g — — + g {E _ Ea)[E _ mE _ Ei) +0{9). (2.9) 

We denote here the states from the ry-subspace by latin letters and the ones from the A-subspace by 
greek letters. Using the relation 

{E - E n ){n\4>) = 0(g 2 )(n\4>) , (2.10) 

we can substitute E in the denominator of the last term in eq. (2.9) by E n . Thus we can now express 
cq. (2.7) in the form: 

(» 2 %ir + * vS?-g)ffU) ) m + ° (/) - (E - Em] {mW ■ (2 - u) 

The first term in the left-hand side of this equation can be expressed in the form 

g2 Hl a Hl n {n ^ = ( 2 Hl a Hl n + g2 Hi Hl n _ E) \ ^ (2 12) 

Hi — h, a \ h, n — £j a (H, — £j a )(£j n — h, a ) J 

Performing an iteration of eq. (2.11) one expresses the last term in the right-hand side of eq. (2.12) as 
a fourth order term and replace E by E n . The final result for the energy independent Tamm-Dancoff 
potential, which differs from the original one in eq. (2.8) by higher order terms, is then 

(V cS ) mn = g * H ]™ H ^ (2.13) 



E n — E, 



a 



ijl ul zrl fjl ul ul ul ul 

4 11 ma n a3 n B~i n ~fn 4 n ma n al n H n 7n 

+ 9 77? 7^T7^ 7^T7^ JTT — 9 



(E n — E a )(E n — Ep){E n — E 7 ) (E n - E a )(E[ — E a ){E n — E y ) 

Even after elimination of the explicit energy dependence, eq. (2.7) does not correspond to the ordinary 
Schrodinger equation since the resulting effective Hamiltonian turns out to be non-hermitean. The 
problem is that the projected nucleon states \<p) have a normalization different from the the states 
we have started from, which are assumed to span a complete and orthonormal set in the whole Fock 
space: 

(frlh) = (vEM^) - MV,-) = 5ij - (frlHjX ( E _\ HA ) ^Hih) ■ (2.14) 

Note that the components ipi in this equation do, in general, not vanish. Thus in order to end up 
with the ordinary Schrodinger equation one has to switch from the set {(/)} of eigenstates of eq. (2.7) 
which is assumed to be complete in the jy-subspace* 9 but is not orthonormal, to another complete 
and orthonormal set {x}- This has been achieved in an elegant way by Fukuda et al. [20] and 
Okubo [21]. They introduced an operator A = XAr] which maps the states 4> into the states 

m=A\4>), (2.15) 

so that the state $ is given by 

m = (l + A)\<p). (2.16) 



* 9 This is presumably true at least in the case of weakly interacting fields. 
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The operator A is equivalent to the operator which enters eq. (2.6) with the only difference of being 
energy independent. One can show that it has to satisfy the equation [21]: 



X(H - [A, H] - AHA) 7] = . (2.17) 
It is now obvious that the states |%) 

| X ) = (1 + ^)1/2^) (2 . 18) 

lie in the r/-subspace and are normalized properly. The effective Schrodinger equation for \x) can now 
be derived from eqs. (2.5), (2.15), (2.17) and (2.18), see ref. [21]: 

(H + V cS )\x)=E\x), (2.19) 

with the effective potential given by 

V cS = (1 + A^AY l l 2 (r, + A^)H{r, + A)(l + A)' 1 ' 2 - H . (2.20) 

A more elegant way to end up with the same effective Schrodinger equation has been proposed by 
Okubo by using a unitary transformation of the form 

/ ^1+1^-1/2 -It(l + IIt)-l/2 \ 
U ~ { I(l+iti)-V2 A(l + iit)"l/2 J, {^) 

with the operator A = XArj. The operator A has to satisfy the decoupling equation (2.17) in order 
for the transformed Hamiltonian to be of block-diagonal form and can (under certain circumstances) 
be identified with the operator A from eq. (2.15). In [2] we have shown how to solve the decoupling 
equation (2.17) and to derive the nuclear force according eq. (2.20) using the methods of CHPT. In 
what follows we will perform a complete renormalization of the NN potential at NLO and demonstrate 
that the renormalized expressions for the potential in this formalism (with the specific choice of the 
unitary operator given in eq. (2.21)) agrees with the one obtained by using the S-matrix methods. 
We will also show how to perform renormalization in the method of unitary transformation and how 
to include pion tadpole contributions and recover the same expressions for renormalized quantities 
as found in covariant perturbation theory using the technique of Feynman diagrams. It is important 
to perform a complete renormalization of the effective Hamiltonian and to take into account all the 
implicit dependence on the quark (or, equivalently pion) mass. 

Before closing this section we would like to stress that the effective Hamiltonian derived with the 
method of unitary transformation is by no means a unique one, as it has also been the case in the 
S-matrix method. The unitary operator (2.21) is not the most general one. Indeed, it is always 
possible to perform additional unitary transformations in r] and A-subspaces, which are not taken into 
account in the definition (2.21). The above unitary transformation is in that sense the "minimal" one. 
This ambiguity with respect to the effective Hamiltonian clearly does not mean that the physically 
observed quantities are defined ambiguously. In what follows, will not be concerned any more with 
these subtleties. 



2.2 One-pion exchange contribution to the scattering amplitude 

The contribution from the OPE to the NN force shown in fig. 1 appears (in Weinberg's power counting) 
at leading order (LO) in the chiral expansion together with two different contact interactions without 
derivatives. The OPEP gets renormalized at NLO due to pion loops and counter term insertions. As 
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Figure 1: Leading order one-pion exchange. The solid (dashed) lines refer to nu- 
cleons (pions). The heavy dots are leading order vertices from C^ N . The diagrams 
resulting from interchange of the nucleon lines are not shown. 



discussed in the previous section, the renormalized OPEP in the S-matrix method is defined by off— 
the-energy shell extension of the field-theoretical amplitude. The corresponding diagrams are shown 
in fig. 2. The graphs 1-4 lead to the renormalization of the nucleon lines, while the diagrams 10 and 
11 contribute to renormalization of the pion line and 7-9 renormalize the pion-nucleon coupling. The 
contribution from graphs 5 and 6 involves an odd power of the loop momentum I to be integrated over 
and thus vanishes. The underlying chiral Lagrangians are given by: 



£(2) 
£(4) 



'-■kN 



,(2) 
-ttN 



c 



(3) 
■kN 



^(V^V M [/t + x+ ), 

^3 , V 2 

Tq (x + ) 



(2.22) 



+ ^ { 2(^*7^) <x+) + 2( X ^U + X UhU ] ) - 4(x t x) - (x+) 2 } 



+ 



N„ 



N„ 



i(v-D)+g A (S-u) 



N„ 



1 



1 



-D.D + Cl ( X+ ) 



N v + 

.2 m 2m 
d w S ■ u{ X+ ) + id 18 S^[D^ x-} + d 28 (i{x+)v ■ D + h.c.) 



N v + 



where F, Q a and m refer to the bare pion decay constant, the bare nucleon axial vector coupling 
constant and the bare nucleon mass* 10 and c\ and d. L are low-energy constants. The brackets ( ) 
denote traces in the flavor space. We adopted here the heavy baryon formulation for nucleon fields, 
with and = (l/2)i^a liu v u denoting the four-velocity and spin operator, respectively. The 
unitary 2x2 matrix U in the flavor space collects the pion fields and is defined in the cr-model gauge, 
which we will use throughout this work, as: 

U = 4 \VF 2 - TV 2 + ir-iv . 



F L 



Further, 



u 



i(u* d^u 



(2.23) 



(2.24) 



-[u^d^u], x± = U ^X^ ± ux'u , x = 2B(s + ip), 



* 10 If dimensional regularization is used, F, Sa and m coincide with the corresponding constants in the chiral limit. 
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3 



4 



10 



11 



Figure 2: Diagrams contributing to the renormalization of the OPE. The heavy dots 
are leading order vertices from and C^l while the solid rectangles correspond to 
vertices from C^ N (in case of diagrams 3 and 4 also from C^ N ). The solid diamond 
represents insertions from For remaining notation see fig. 1. 



where the quantity B is just a constant. We do not take into account external vector and axial-vector 
fields. The physically relevant values for scalar and pseudoscalar sources s and p are 

s(x)=M, p(x) = 0, (2.25) 

where M denotes the quark mass matrix. Note that we have only shown explicitly the terms in 
eq. (2.22) we will use in our calculation. We also note that 1/ m terms are suppressed in Weinberg's 
approach and would appear one order higher. Further details on construction and structure of the 
effective Lagrangian can be found in [22, 23, 24, 25, 26]. 

The OPE at leading order corresponds to the diagram shown in fig. 1 and is given by 

9 A 2 1 
A OP e = ^ ( Tl • r 2 ) ^— — 2 (S, ■ q)(S 2 ■ q) . (2.26) 

Here q denotes the momentum transfer of the nucleon, i.e. q = p' — p, where p' and p are final and 
initial nucleon momenta. Further, M is the leading term in the quark mass expansion of the pion mass. 
The more convenient expression in the rest-frame system of the nucleons* 11 with = (1, 0, 0, 0) and 

* n We are only interested in the specific kinematic, in which both nucleons move with the same velocity « M and the 
relative momentum is small. 



= (0, 1/2 a), where the cij are the Pauli spin matrices reads 



9a 



OPE 



(2F) 2 q 2 + M 2 



(ai-q) (a 2 -q) + 0(m- 



(2.27) 



We have introduced here the superscript NC (= "noncovariant" ) in order to distinguish between the 
noncovariant notation of eq. (2.27) and covariant one in eq. (2.26). The noncovariant notation is 
appropriate for calculation of the S-matrix based upon the effective Hamilton operator. Note that 
the physical values gA, F w and M w of the axial nucleon coupling, pion decay coupling and pion mass 
can be used in the expression for the OPE if one restricts oneself to the LO analysis. 
We will now evaluate the NLO corrections to the OPE, which result from the diagrams shown in fig. 2. 
The calculation can easily be performed using e.g. the Feynman rules from [22, 25]. The diagrams 
1,...,4 renormalize external the nucleon lines and lead to 



"4-1,2,3,4 = ((Zn) 2 - l) -A.QPE ■ 



Here Zn is the nucleon Z-factor given by 



1 + S'(0) . 



(2.28) 



(2.29) 



The nucleon self-energy E(tD) (to = v • k, /c-nucleon momentum) is given up to order Q 2 (Q refers to 
a generic small momentum scale) by 



E(£>) 



AM' 



(ci + 2d 28 u) + J 



dH 



(2vr) 4 v {k-l) + iel 2 - M 2 + ie 



(S-l)(-S-l). (2.30) 



We will now proceed in two ways. First, we apply dimensional regularization and the same subtraction 
as in ref. [24] and recover the standard results for the amplitude. Secondly, in order to compare our 
results with the ones obtained in the method of unitary transformation we will go into the rest-frame 
system with = (1, 0, 0, 0) and = (0, 1/2 a), where the cjj are the Pauli spin matrices and perform 
integrations over zeroth component of the loop momentum I via contour methods. Furthermore, we 
will not restrict ourselves to any specific regularization scheme in that case since we are only interested 
in the general structure of the amplitude. The divergent integrals are always to be understood as being 
regularized by some standard methods. In order to distinguish the results obtained in this way from 
the ones found using dimensional regularization and the subtraction defined in ref. [24] we will denote 
the latter ones by the superscript DR. 
For the nucleon Z-factor Zjy in eq. (2.29) we find: 



7 DR 



where we have introduced 



3 9 A 2 M 2 
32tt 2 F 2 



M 

3 In— + 1 

A 



I -3^J 13 -8M 2 d 28 . 



8F 2 



J n 



d 3 l I 



2m 



(2vr) 3 UJ r , 1 ' 



-8M 2 d r 28 (X), 



(2.31) 



(2.32) 



Clearly, Z N R can also be obtained directly from Zn performing dimensional regularization in three 
dimensions of the divergent integral J13 . 
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The nucleoli mass shift is related to the self energy via 



ton =m +5mN , 5mN = S(0) . (2.33) 
Using dimensional regularization one finds 

^ = -4 Cl M 2 - , (2.34) 

while the unregularized expression (after performing integration over the zeroth component Iq of the 
loop momentum) reads: 

Sm N = -4 Cl M 2 - ||f J 12 . (2.35) 

Notice that the mass of the nucleon does not appear explicitly in the expressions for the NN potential 
at NLO but enters the two-nucleon propagator in the Lippmann-Schwinger equation. We will discuss 
this issue later on. 

For diagrams 7,8,9 in fig. 2 we obtain in dimensional regularization (D.R.) 

"4^9 = -^w>A)pe, (2.36) 



dr 2M 2 L(X) , M 



^7 8 — ^ r 



7,8,9 ' p 2 T l 67r 2 F 2 



,° 2 „m M ° 2 

{9 A 2 + 2) In — + 5a 2 



+ 



8M 2 <q 6 (A) _ 4M 2 di 8 

o o 

9a 9 A 



where we used the same definitions for L(X) and d^ 6 (A) as in [24]. Note that the constant d±s, which 
leads to the Goldberger-Treiman discrepancy, is finite and does not get renormalized. Switching 
to the noncovariant notation and performing integration over Iq we find the following unregularized 
expressions for these quantities: 

•4-7,8,9 = ^7,8,9 -4oPE > (2.37) 

y , 2 \ 1 / , 8M 2 d m 4M 2 d 18 

^7,8,9 - -^{i-d+J^l) + + Tj— . 

where d = 4 in the physically relevant case of (3+1) space-time dimensions. 
Finally, renormalization of the pion line in graphs 10 and 11 in fig. 2 leads to 

A 10 ,n = [SZ n + 2 _ M 2 ) Ajpe , (2.38) 

where 

2M 2 

Z w = l + 5Z w , 5Z w = --j^ fT^i (2.39) 

Ml = M 2 + 5M 2 , 5M 2 = M 2 {^ + ^ 3 }, 

Here, M w is the physical value of the pion mass. The last term in the brackets in eq. (2.38) gives the 
correction to the one-pion exchange due to the shift in the pion mass and should be ignored if the 
physical value of the pion mass is used in the expression (2.26) for the OPE at LO. The quantity A^ 
is given in D.R. by 

A° R = 2M 2 (l(X) + JL In + 0(d - 4) , (2.40) 
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and in the previously introduced notation with divergent integrals Jjj just by A w = 1/2 J 01 . Note 
that in the calculation of the pion wave function renormalization and mass shift one has to take into 
account the noncovariant pion propagator of the form: 

ci 4 x(0|T^7r a (x)^7r 6 (0)|0) = ^1^-— - ig^g vQ , (2.41) 

as discussed in [29]. Such noncovariant pieces in the pion propagators (as well as noncovariant pieces 
in the interaction Hamiltonian arising from time derivatives) are usually omitted if one works in 
dimensional regularization. This does not lead to wrong results since the missing terms correspond 
to power-law divergences, which vanish in D.R.. However, in a general case such noncovariant pieces 
should be taken into account. For instance, if the last piece in the above expression is omitted and 
cut-off regularization is used, pions become massive in the chiral limit due to the tadpole diagram. 
Summing up the contributions of all graphs in fig. 2 we obtain the renormalized OPE as 

,ren 1 ( 9aV ( ' <iM%d 18 \ T X ■ T 2 _ . . 

where gA/F-n is given in D.R. by 

9A _ 9 A ( g\Ml M n g\Ml 2 ~ AM$ ^ A ^ 

F^-Tl 1 - i^Ff ln X " 16^Fl ~ 8M ^28(A) + —dwW ~ j^UW ) (2.43) 

o 

Note that we substituted 9a and F w in the brackets in the right-hand side of above equations by 
their renormalized values qa and F w , which is correct at this order. If one is only interested in the 
real world observables, which correspond to the physical value of the pion mass, one can equally well 
rewrite this expression in terms of scale independent LECs d\Q and I4 as follows: 

9a_°9a( a g\Ml AMI Ml \ 

The LECs die and I4 are defined according to refs. [24] and [6]: 

^e(A) - dl6 -^L(A) = J 16 -^lnA ( 2 .45) 

ll(X) ee k-2L(X) = T ^h + ^ln^, 

where k±q =gA (4— gA )/8. The renormalized LEC cfos has no finite piece (or, in other words, the 
corresponding scale-independent LEC J28 vanishes): 

4(A)^ 2 s-ff^A) = -^lnA, ( 2 .4 6 ) 

with K28 = —9/16 g°A 2 ■ Thus the corresponding counter term in the first line of eq. (2.31) cancels 
against the logarithmic term. We note that it is a convention that the contact terms needed for 
renormalization do not have a finite part. In principle, one could allow for such contributions but that 
would only lead to a reshuffling of the finite parts of other LECs. If one insists on working with the 
minimal number of independent terms, such a choice has to be taken (as discussed in more detail in 
[24]). 

Finally, we display the unregularized expression for the ratio gA/F w : 

9a_°9a(, g\( A , , 2\ M2 su2L , 4M, 2 



F-k F 



(l - (4 - d + ^) J U - - SM% 8 + ^) . (2.47) 
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2.3 One— pion exchange from the method of unitary transformations 

Let us now calculate the renormalized OPEP using the method of unitary transformation, following 
the lines of refs. [2, 3]. The essential difference compared to our previous work [2] is that we now have to 
take into account pion tadpole contributions. The corresponding diagrams lead to renormalization of 
various parameters in the effective Lagrangian and were of no importance in the previous calculations 
[3], where physical values for renormalized coupling constants have been used. The pion tadpole 
contributions are, however, important for the present analysis, since they lead to an additional pion 
mass dependence of various renormalized quantities. In the following we will show how to deal with the 
pion tadpoles and how to perform complete and consistent renormalization in the method of unitary 
transformation. We would also like to stress that renormalization of the effective Hamilton operator 
within the method of unitary transformation has been considered recently in a different context, see 
refs. [27, 28]. 

Let us begin with the effective Lagrangian C wn for pions given in the first two lines of eq. (2.22). The 
relevant terms in the cr-model gauge, eq. (2.23), read: 

4 2 2 = \ (d,7T ■ &>* - MV) + JL (^TT • 7T) 2 - ^M 2 7T 4 + . . . (2.48) 

We now introduce the renormalized pion field and mass: 

vr r = Z~ 1 ' 2 TV, Z 7T = 1 + 5Z 7T , (2.49) 
Ml = M 2 + SM 2 , 

where SZ n , (SM 2 )/M 2 ~ C(Q 2 /A 2 ). We express eq. (2.48) in terms of the renormalized pion field 
and mass by simply replacing tv — > 7r r , M 2 — > M 2 and 

t^kn > cSve ~\~ S (^nl^j where 

5 (4 4 j) = Uz^ r ■ d^ r - \Ml (bZ^ - 7T 2 + . . . . (2.50) 

Here the ellipses refer to terms with more pion fields. In what follows we will always work with 
renormalized pion fields and therefore omit the subscript r. 

The Hamilton density corresponding to the Lagrangian introduced above can be derived using the 
canonical formalism and following the lines of [29]. Let us first express the effective Lagrangian in a 
more compact form: 

1 M 2 

C*n = -d^ a G ab {7T)d^ b - -flT 2 B(ir) , (2.51) 

where 

(11 M 2 \ 1 
l + 5Z n + — ^ J + -j^iT a ir h + ... =5 ab + G ab (n), (2.52) 

B(rt 1 4 17 4 2kM " 4 2UM " 4 1 it 2 4 14. R(tt) 

B(7V) - 1 + SZ V - -j^- + — + — + ^2 7T +... - 1 + B(7V). 

Here the subscripts a, b correspond to isospin indices. For the momentum n a conjugate to the field 
ir a we get: 

U a = ^ = G ab (7r)d 7r b . (2.53) 
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The Hamiltonian H n7T is given by 



1. 



Ml 



= Ii a G-J(-K)Ii b - C m = -U a G- b \7v)n b - -dutaGriW&irb + -f^B(v) 
We now divide the Hamiltonian into its free and interaction part as follows: 



(2.54) 



7T7T ' 

' 1 M 2 

H° n7T = -U a U a --d i 7T a d i 7r a + ^7r 2 



(2.55) 



n 



-\d^ a G ab (Tr)d^ b - ±d ir a G 2 ab (ir)d ir b + ^M 2 tv 2 B(tv) . 



Since the field operators 7r a , <9j7r a and Il a transform simply into the interaction picture (IP), ir a — > 7r* p , 
di7T a — ► (5j7r a ) IP and Il a — ► (<9o7r a ) IP , we can easily express the interaction Hamiltonian 7^ in the IP 
representation as:# 12 

K n = -i^7r a G a6 (7r)^7r 6 + iM> 2 B(7r) + ia 7r a {G' 2 (7r)(l + G(7r))- 1 ) d n b (2.56) 

Z Z Z I ) ab 

5M 2 t 2l 3 M 2 2l A M 2 



(7r.a At 7r)(7r-^7r) + -^M 2 7r 4 + ..., 



F 2 



F 2 



2F 2 y " M v v " " ' 8F 2 ~ 
where we have omitted the superscripts IP. We now write TT a (x) in terms of creation and annihilation 
operators a a (k), a a (k) as 

d 3 k 1 



J (2tt) 6 i z yJ2u k L 

where io^ = J k 2 + M 2 and the operators a Q (fe) and ai(&) satisfy the commutation relations 



a a (fc), a 6 (fc') = al(k), a\{k') 



0, 



a (fc), aj(fc') = <5 ab o(fc-P). 



(2.57) 



(2.58) 



Substituting eq. (2.57) into eq. (2.56) and performing normal ordering, we end up with terms of the 
form a^a, aW , aa, a)a)aa, .... In fig. 3 we show symbolically the contributions to the Hamiltonian 
of the form a^a, a^a^ and aa. Note that the closed loops result from contraction of the operators a 
and a' in the terms in the last line of eq. (2.56) when bringing these into normal ordered form. For 
example, for the term of the form a^a we get the following expression: 

HL = £ / d * k ^ -lCk)a a (k){h 2 (-5Z„ - *0 - JL 



(2.59) 



+ ( <^ 



M^ 2 
M 2 



2/ 3 M 2 2/4M, 



+ 



F 2 



+ 



F 2 



+ 



4F 2 



■^01 



+ ... . 



We now require that there are no terms of the form a^a, aW and aa in the interaction Hamiltonian at 
the order considered in the chiral expansion (since we already work with the renormalized pion fields) . 
Then we obtain from eqs. (2. 49), (2. 59): 



Z* = 1- 



2hMl 
F 2 



2F 2 



Joi , 



(2.60) 



Mi 



M 



1 + ^ + ^ Jn 



F 2 



4F 2 



* 12 In the last line of eq. (2.55) one has first to replace <9o7r a by G a ^(Tv)Ii b 
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Figure 3: Contributions to the one-pion Hamiltonian. The heavy dots are leading 
rder verl 



order vertices from C^l while the solid rectangles correspond to vertices from d£l + 



These expressions coincide precisely with the ones found using the S-matrix method, eq. (2.39). 
Let us now include nucleons. The relevant terms in the effective Lagrangian C n N from eq. (2.22) are 

4n = A ft ^ -|^ra.V7r-^(7r ( T-V7r)r.7r)^iV + ..., 

4n = Nl (sm N + + 4 Cl M^j N + ... , (2.61) 

£ 7riv = N I I p — + — p Ip" 5 ^ I tct • Vtt + 8id 28 M n d J N + ... . 

Here several comments are in order. First, the superscript (i) in the above Lagrangians refers, as in 
eq. (2.22), just to the number of derivatives and/or pion mass insertion. Secondly, we prefer to use 
the physical and not the bare mass of the nucleon in the heavy baryon expansion, i.e. we parameterize 
the four-momentum of the nucleon p^ as p^ = m^v^ + Z^, where is a small residual momentum 
v ■ I <C itin- Expressing the Lagrangian in terms of the physical instead of the bare nucleon mass is, of 
course, just a matter of convenience and all results remain unchanged. The physical and bare nucleon 
masses are related by 

rriN =m +5rriN , dm^ = Smff + Smff + . . . , 6m$ ~ O (Q l ) . (2.62) 

Note that 5mff = —Ac\M^ as can be seen from the second line of eq. (2.61).* 13 As in the previous 
section we restrict ourselves to the rest-frame system of the nucleons with = (1,0,0,0). Notice 
that the term N^5m^N enters the Lagrangian C^ N as a consequence of using the physical nucleon 
mass instead of the bare one. Further, we have not included the term — l/(4F 2 )A^r • (n x -k)N in the 
effective Lagrangian 4jv'' wri ich results from the covariant derivative of the nucleon field in the third 
line of eq. (2.22). As already explained above, the corresponding contributions to the OPE vanish since 
an odd power of the loop momentum / to be integrated over enters the explicit expressions. Finally, 

* 13 Pion loops start to contribute at next higher order in the chiral expansion. 
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it should be kept in mind that the effective Lagrangian (2.61) is expressed in terms of renormalized 
pion fields, which leads to the additional contribution with an insertion of 5Z^ in C^ N - The Hamilton 
density corresponding to eq. (2.61) is given by 

V 2 

n\ = -n ] - — n , 



2m n 



Ho = N f ( |^TCT-V7r + ^(7r«T- V7t)t-7t j N + ... , (2.63) 



to ■ VtvN + ... . 



Here we switched to the field N defined as 



N = (Z^-^N =fl- -JLj = (1 + AM^d 28 )N , (2.64) 

in order to get rid of the time derivative in the third line of eq. (2.61). Further we have used a 
classification of the vertices in effective Lagrangian different from ref.[2]. To be more precise, the 
dimension A of the interaction is defined as: 

A = d+^n + l-2. (2.65) 

Here d, n and I is the number of derivatives or insertions of M n , nucleon field operators and insertions 
of the factors l/m^v, respectively. According to the power counting suggested by Weinberg [1] we treat 
the nucleon mass as a much larger scale compared to A x (m^ ~ A 2 /Q). As a consequence, relativistic 
corrections enter at higher orders compared to the corresponding chiral corrections. Notice further that 
the nucleon field operators in eq. (2.63) are always taken in the normal ordering. Allowing contractions 
of nucleons field operators would only lead to shifts in the values of various coupling constants in the 
effective Lagrangian (Hamiltonian) which are independent of M T and thus of no importance for our 
considerations.* 14 

We will now follow the lines of [2] and derive the effective Hamiltonian, which acts on the purely 
nucleonic subspace of the whole Fock space, using the method of unitary transformation. The idea of 
this method is described in section 2.1. For a detailed discussion on the way of solving the nonlinear 
decoupling equation (2.17) and calculating the effective Hamiltonian according to eqs. (2. 19), (2. 20) 
the reader is addressed to ref. [2]. 

Let us begin with the one-nucleon effective Hamiltonian -Sin: 

H m = + + . . . , (2.66) 



MO) 



ill ! H 2 - H ^—Hq ) in 



* 14 Remember that one does not start from a normal ordering for pion field operators is required. Pion tadpole diagrams 
lead to M^-dependent renormalization of various constants in the effective Lagrangian and thus must be considered 
explicitly. 
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where the subscript i (superscript j) of the projector r\i (A- 7 ) denotes the number of nucleons (pions) 
in the corresponding state and u> is the pionic free energy. Further, Hi denote Hamilton operators 
corresponding to the Hamilton density Hi in eq. (2.63). The operator X 1 A r]i is given by [2]: 

X^om = Hovi ■ (2.67) 

The superscript v of refers to the order in the chiral expansion and is defined as:# 15 

v = -2 + 2E n + 2(L ~C)+Y^ ViAi , (2.68) 

i 

where E n , L, and C denote the number of nucleons, closed loops and separately connected pieces, in 
order. Furthermore, Vi is the number of vertices of type i. Performing straightforward calculations 
we end up with the following result: 



- P 2 Xm^ 3 9a T (0 fiQi 

h ™-2^~ 6mN ~ 8i^ Jl2 ' (2 ' 69) 



where we switched to a more convenient nonrelativistic notation, i.e. we omite creation and annihila- 
tion operators as well as summation over the appropriate quantum numbers. To be more precise, the 
relation between if in and hm is given by 

H m = I d3 P n *(P) MP) hm{p) , (2-70) 

a J 

where a refers to discrete quantum numbers (spin and isospin) and ria(p) (n a (p)) is a creation 
(annihilation) operator for the nucleon field. It is now easy to read off the expression for the nucleon 
mass shift Sm^ up to the order considered here: 

5m N = -4 Cl M2 - ^|f-J 12 , (2.71) 

which agrees precisely with the result obtained in the previous section using the covariant perturbation 
theory. One sees from eq. (2.69) that the leading one-nucleon Hamiltonian is given, in accordance 
with our expectation, just by the nucleon kinetic energy = p 2 /(2mjv). Notice that since we are 

(2) 

interested in the NN interaction at order u = 2, we should in principle also consider and thus 
include higher order nucleon mass shifts. We will, however, see explicitly in what follows that this is 
not necessary. Moreover, we will show that even the leading shift in nucleon mass, 
contributes at NNLO and thus can be neglected. 

The nucleon Z-factor Zn needs, in principle, not be calculated separately in the method of uni- 
tary transformation, since it already enters eq. (2.20). More precisely, the part of the Z-factor 
ZJj, which corresponds to dressing of the bare nucleon by virtual pions, is related to the operator 
771(1 + A)~ 1 / 2 r]i. Indeed, the Z-factor is given by 

\ Xl ) = (Z^)- 1 / 2 ^) , (2.72) 



# 15 We use here the definition of the counting index v different from the one introduced in ref. [2]. The presently used 
notation is more transparent when operators with different number of nucleons, i.e. f N, 2N, 3N, . . . forces, are considered. 
All results of ref. [2] remain however unchanged. 
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where \4>i) is the original (bare) one-nucleon state and |xi) is the properly normalized one-nucleon 
state (i.e. which satisfies the orthonormality condition as described in section 2.1). It follows from 
eqs. (2. 72), (2. 18) that: 

{Zl)- 1 ' 2 = (0 1 | Xl ) = ^|r ?1 (l + AU)V2 r?1 | 0) (273) 
= (<f,\ m (l + \a^Aq + . . .) V i\4>) = 1 + ^Jis + ■■■■ 

The complete nucleon Z-factor Zn, which is related to the specific choice of nucleon field made in the 
Lagrangian (2.61), is given by (at the order considered) 

2 

Z N = Z%Zj r = l-8M%dv i -^gj 13 . (2.74) 

This agrees with Zn from eq. (2.31). We should, however, stress that the Z-factor is not an observable 
quantity and depends on the choice of fields. For example, there would be no contribution proportional 
to d,28 in Zn if we had decided to work with the fields N instead of N. For more discussion on the 
role of the nucleon Z-factor in the old-fashioned perturbation theory and Hamiltonian formalism the 
reader is referred to [21, 20, 31, 32, 33] (for corresponding discussions in the framework of baryon 
CHPT, see e.g. [34, 35]). 

Let us now consider the 2N effective Hamiltonian at order v = 2. The explicit operators at NLO are 



Figure 4: Contributions to the two-nucleon Hamiltonian at order v = 2: discon- 
nected graphs. Heavy dots refer to the leading order ttNN vertex of dimension 
A = while the solid rectangles correspond to insertions of the vertex with A = 2. 
For remaining notation see fig. 1. 



given in ref. [2]. For our purposes we do not need all terms but only the following ones: 

^2N = — Ho Hq ; Hq Hq + -Hq——^ HqT] 2 Hq Hq + -Hq Hq^Hq-^-^Hq 

\ UJ U)\ + UJ 2 W 2 (LO) z UJ 2 LO (LO) z 

+ Hq—H 2 —H - ±Ho-^-Hq V2 H 2 - }-H 2m Ho-^-H ) m . (2.75) 
uj uj 2 [ujy 2 (uj) z J 

In fig. 4 we show disconnected diagrams, which contribute to the 2N effective Hamiltonian at order 
v = 2, which do not include the nucleon self-energy.* 16 Denoting by Ai the common spin-isospin 

* 16 Notice that the third diagram in fig. 4 is related to the last two operators in the first line of eq. (2.75) and thus does 
not correspond to the iteration of the nucleon self-energy contribution. The same holds true for last two graphs in this 
figure. 
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structure as well as the loop integrals and using the first line of eq. (2.75) we see that the contribution 
of the first three diagrams in fig. 4 vanishes: 

/ 2 1 1 _L_ , 1 

\ u>fu>2(ui + ^2) w^o;2(a;i + ^2) wicjf^i + ^2) ^1^2 ^1^2 

A similar cancellation is also observed for the remaining graphs in fig. 4 using the second line of 
eq. (2.75). Thus no disconnected diagrams (which are different from the nucleon self energy graphs) 
contribute to the NN potential at NLO. 

Consider now the diagrams shown in fig. 5, which are of the same topology as Feynman graphs 1-4 in 
fig. 2. The corresponding operators are given in eq. (2.75). Note that the contribution of the diagrams 



J = (2.76) 



Figure 5: Contributions to the OPE at order v = 2. For notation see fig. 4. 



in the lower row in fig. 5 vanishes similarly to the previously considered case of disconnected graphs. 
Performing explicit calculations one finds the same result as in the covariant approach, namely:* 17 

h^ = Z N A^ E , (2.77) 

where Zjy is given in eq. (2.31). The OPE amplitude A^ E is the same as -4qpe ln (2-27) with 
the only difference that the physical pion mass M n is used instead of M (since we are working now 
with renormalized pion fields). We will not consider the topologies related to graphs 5 and 6 in fig. 2, 
since the corresponding contributions vanish as explained above. The remaining contributions to the 
OPE at NLO in the method of unitary transformation are depicted in fig. 6. The first four diagrams 
show all possible time orderings of the covariant graph 7 of fig. 2 and correspond to the first operator 
in eq. (2.75). Explicit evaluation leads to the same expression as found in the covariant approach, see 
eq. (2.37): 

h^ = Z 7 A^ E , Z 7 = -|iJ(l-d+^y J 13 . (2.78) 

The two pion tadpole diagrams in fig. 6 result from contracting two pion field operators when expressing 

o 

the term —N^ -^(na ■ V7r)r • irN from eq. (2.63) in normal ordered form, as it also happened in case 

* 17 In order to be consistent with the results of the previous section we give here and in what follows explicit expressions 
for the effective Hamiltonian based upon the nucleon field N, which is related to N via eq. (2.64). As a consequence, we 
have to take into account the additional normalization factor (Z^)^ 1 ^ 2 . 
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Figure 6: Remaining contributions to the OPE at order v = 2. For notation see 
%• 4. 



of the pion self-energy discussed above. A straightforward but somewhat tedious calculation leads 
again to the same result as in the covariant approach: 

^2N = ^8-4oPE > Z S = ^01 • (2.79) 

Consider now the last two diagrams in fig. 6. Since we have chosen to work with renormalized pion 
field operators in the method of unitary transformation, we have no diagrams corresponding to graphs 
10 and 11 in fig. 2. Contributions of these two graphs are now taken into account by using the physical 
pion mass M n instead of M as well as by an additional interaction proportional to SZ W in the last line 
of eq. (2.63). For the last two diagrams in fig. 6 we thus find: 

4*>=(^ + ^-^W E . (2.80) 

Summing up all corrections to the OPE in the method of unitary transformation we end up with the 
same expression (2.42) as found from covariant perturbation theory using the S-matrix methods. 
In summary, we have shown how to perform renormalization in the method of unitary transformation. 
We recover the same expressions for pion and nucleon mass shifts as well as for corresponding Z- 
factors. We have further demonstrated that the NLO OPE potential derived using the method of 
unitary transformation, where the unitary operator is chosen according to eq. (2.21), coincides with 
the off-the-energy shell extension of the S-matrix. The TPE contribution has already been considered 
in both approaches, see refs.[16, 2, 3] and identical results have been reported as well. We do not need 
to reconsider the TPE contributions and the pion loop diagrams which renormalize the short-range 
terms since no pion tadpoles appear in those cases and the results of refs.[2, 3] can be adopted without 
any changes. 

2.4 Explicit expressions for the potential at NLO 

We now give the explicit expressions for the chiral effective NN potential VnlO) which we use for 
extrapolation in the pion mass: 

Vn L o = V OPE + V TPE + V coat , (2.81) 

where 

f ope _ JAL + 2A _ ^i s \ Tl . T2 (hJ^Ji _ (M2) 

4 F% \ 9 A q 2 + M% 



20 



yTPE = 



'U%(5g\ - Ag\ - 1) + q\2Zg\ - 10g\ - 1) + 



4M 2 +q 2 } 

+ ^1^(23^-10^-1)} (2.83) 

"eSfe ( L(9) + ln ^) {^■^■^-^^■^}, 



+ dq 2 + C 2 fc 2 + {C 3 q 2 + C 4 fc 2 )(ai • ct 2 ) 

+ iC 5 ai ^ • (fe x g) + C 6 (<f ■ B\){q- a 2 ) + C 7 (£ • ai)(k ■ a 2 ) , 

with g\ and the physical values of the nucleon axial coupling and pion decay constant, respectively. 
Here and in what follows we denote the value of the pion mass by M v in order to distinguish it from 
the physical one denoted by M w . Further, 



'AMI + q 2 \UM 2 + q 2 + \q I 
L(q) ee L(\q\) = ^— ln V , (2.85) 

i y i Ziivi'fl 
and A represents the relative shift in the ratio gA/F^ compared to its physical value: 

A = (W^)m, ~ (9A/F n ) M ^ 

{9 A /Fir) M n ' 



Note that in the TPEP we only take into account the explicit M^-dependence and use the physical 
values for g& and F n . This is perfectly sufficient at NLO since any shift in g& and F n for a different 
value of in the TPE is a N 4 LO effect. We have expressed the potential in eqs. (2.82)-(2.84) in 
such a way that it coincides for M n = M n with the one given in [3].# 18 The constants Cs,t and Ds,t 
are related to the Cs,t from [3] via 

C s ,t = C s ,t + M 2 D s ,t ■ (2.88) 

Note further that the short-range terms of the type M 2 lnM^ in eq. (2.84) result from the two-pion 
exchange as well as from the renormalization of the leading-order contact forces by pion loops. It 
is important to stress that renormalization of the LECs C5, Ct, Ci v ..7 due to pion loops does not 
depend on the pion mass and thus is of no relevance for this work. 

We will now briefly discuss an important issue related to the renormalization of the NN potential. We 
have performed the calculation of the renormalized OPE in the previous section without specifying 
a regularization scheme and gave only general (unregularized) expressions. One might worry about 

* 18 The term in eq. (2.82) which breaks the Goldberger-Treiman relation and is proportional to dis is not shown explicitly 
in ref . [3] . 
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the dependence of the results on the regularization scheme. Let us consider as a simple example the 
nucleon mass shift in eq. (2.71). Applying dimensional regularization to the divergent integral Ju one 
ends up with the finite shift given by eq. (2.34). The bare nucleon mass m in that case is finite and 

o 

the constant c\ is not renormalized. Similarly, some other bare quantities like M, F, 9a are finite. 
The situation might change if a different regularization scheme is applied. For example, performing a 
momentum cut-off regularization of Jyi in eq. (2.71) one finds: 

Sm N = -4 Cl M* - ||f J 12 = -A Cl Ml - ||f (aA* + 0M*A + ^ + ©(A" 1 )) , (2.89) 

where A is the momentum cut-off and the coefficients a and (3 depend on the precise choice of the 
regulator. Taking the limit A — > oo and performing a redefinition of the constants m, c\ in order 

o2 

to absorb the terms proportional to A 3 and A, we obtain the same finite shift -3 9 A M 3 /(32vrF 2 ) 
as in the case of dimensional regularization. The difference to the previously discussed case is that 
the bare nucleon mass and the bare constant c±, which do not correspond to observable quantities, 
are now infinite. Also other bare parameters are infinite if the cut-off regularization is used, see e.g. 
refs. [36, 37]. It is, however, crucial to understand that, say, the shift in the nucleon mass from its 
observed value to the one in the chiral limit is a well defined and finite quantity and is not affected 
by a specific choice of the regularization scheme. The same holds true for F n and qa- Therefore, we 
are free to adopt the results for the corresponding shifts in these constants found in CHPT analyses 
of various processes and based upon the dimensional regularization (i.e. we can use the values for the 
LECs in eq. (2.86)). We will specify the numerical values of the appropriate LECs in the next section. 
To close this section let us make a comment on calculating various observables based upon the effective 
potential introduced above. The NN phase shifts are obtained using the T-matrix method. The 
corresponding partial-wave projected Lippmann-Schwinger equation reads: 

„ - ,) + E/ ^,(p'. » ■ 

Notice that the nucleon mass enters the expression for the two-nucleon propagator in the above 
equation. Thus apart from the explicit and implicit M^-dependence of the NN force, which has 
been discussed above, one should, in principle, take into account the M„--dependence of mjy in the 
Lippmann-Schwinger equation. Since the nucleon mass itin is treated in the Weinberg power counting 
formally as a much larger scale compared to A x (i.e. tun ~ A^/Q), the shift in the nucleon mass gives 
the correction to the T-matrix which is suppressed by 3 powers of the small momentum scale and 
thus contribute at NNLO: 

m ^ = m CL + 4c , M 2 + =m CL^ + 4^ + ^ =m CL + ^ j (291) 

where m CL refers to the nucleon mass in the chiral limit (=m if dimensional regularization is applied). 
We have checked the above estimation numerically and found indeed a significantly smaller effect 
in the deuteron binding energy (in the chiral limit) from the nucleon mass shift compared to various 
contributions at NLO. It goes without saying that the potential V(p' , p) is multiplied by the regulating 
functions /r(|p|), /r(|p'|) in order to cut off the large momentum components in the Lippmann- 
Schwinger equation (2.90), which cannot be treated properly in effective field theory. We used here 
the same exponential function /r(\p\) = exp[— p 4 /A 4 ] as in ref. [4] and vary the cut-off A in the range 
from 500 to 600 MeV. 
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3 Results 



3.1 Chiral input and contact terms 

For the OPE contribution in eq. (2.82) both explicit and implicit dependences on the pion mass are 
known, so that a parameter-free extrapolation to values M n away from the observed one is possible. 
In what follows, we use g\ = 1-26, F w = 92.4 MeV. The constant U is fixed from the scalar radius of 
the pion and the corresponding value is I4 = 4.3 [6]. The LEC die has recently been determined from 
the process ttN — > ttttN [38]. In the following we will use the updated value from ref. [39]. To be 
more specific, we average the values obtained in the 3 different fits in the last reference. This leads to 
Ji6 = —1.23 GeV -2 . Note that although all values of die from [39] have the same sign (negative) and 
are of a similar magnitude, the uncertainty in the determination of this constant remains quite large. 
We come back to this at the end of this section. The constant d±s is fixed from the observed value of 
the Goldberger-Treiman discrepancy. Using the values of g n N extracted from irN phase shift analysis 
[40], we get the following value for the LEC d 18 : d 18 = -0.97 GeV~ 2 . Thus we obtain for the ratio 
9ttN /mN, which is nothing but the strength of the OPE and is related to g& and F w via (at the given 
order in the chiral expansion) 

_ £ / _ 2Ag 3 \ _ 



the following value in the chiral limit 

-15.0 GeV- 1 . (3.2) 

m N J CL 

Comparing this value with the physical one of 13.2 GeV -1 , we conclude that the OPE becomes stronger 
in the chiral limit. 

We now turn our attention to the remaining contribution in the NN potential (2.82). The TPE part 
is parameter-free and the values of the LECs Ci v ..j as well as combinations Cs,t + M%Ds,t have 
already been fixed from the fit to low-energy data in the NN S- and P-waves. Unfortunately, the 
LECs As,T, which contribute to S-wave projected potential, are not known at the moment. Ideally, 
they should be fixed from the NNir system in processes like e.g. pion-deuteron scattering. This has 
not yet been done. In order to proceed further we assume natural values for these constants, i.e.: 

Ds > T = F 2A2 ' where a S,T ~ 1 , (3.3) 
" X 

and A x ~ 1 GeV. In our analysis in [43] we have shown that all values of the dimensionless coefficients 
a related to the contact terms lie at NLO in the range —2.1 ... 3.2 for all cut-offs employed. Thus 
making the conservative estimate 

-3.0 < a s , T < 3.0, (3.4) 

we are rather confident to be on the safe side. Note that this corresponds to the following range for 
the partial-wave projected values: 

48vr - 48vr 24vr - 24vr 

7T X T X I" X T X 

It is worth mentioning that the term M%Dx(cri • a?) breaks Wigner symmetry. The value of Dt 
expected from eq. (3.5) leads to a Wigner symmetry breaking effect that is comparable in size with 
the one from the observed nonvanishing (leading order) value of Ct [43] (i.e. M^Dt ~ Ct), so that 
the possibility of exact Wigner symmetry in the chiral limit exists. In that case, Wigner symmetry 
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breaking would be related to the explicit chiral symmetry breaking of QCD. Last but not least, we 
use here the values of the LEC Cs,t, Ci,...,7 from [3].# 19 

Before showing our results for the phase shifts and other quantities, we would like to make a comment 
on the uncertainty of our extrapolation in M n . As already pointed out, the main source of uncertainty is 
related to the unknown values of the LECs Dg : T, which enter the expression for contact interactions at 
NLO, see eq. (2.84). Although the LECs I4, di6,i8 are known from CHPT analyses of various processes 
and thus the complete M^-dependence of the ratio g-wN/mN is fixed at NLO, some uncertainty still 
remains due to the uncertainty in the determination of the abovementioned LECs. While the constant 
I4 is known with relative small error bars, so that the corresponding uncertainty needs not to be 
discussed here, the LECs di6,i8 are not known very precisely. Using the values of g w N extracted from 
three different irN phase shift analyses [41, 40, 42] one gets the values d±s = —1.54 GeV -2 , d\8 = —0.97 
GeV~ 2 and d 18 = -0.84 GeV -2 , respectively. From eq. (2.82) one sees that dig does not contribute 
to the OPE in the chiral limit. Thus at first sight our predictions in the chiral limit seem not to be 
affected by the uncertainty in this LEC. In fact, this is not quite true, since changing the value of dis 
leads to a modification of the strength of the OPE at the physical point (M n = M n ). Therefore, one 
should, in principle, refit at the same time the LECs Cs,t, Ci,...,7 in order to describe the corresponding 
phase shifts. Modifying the LECs related to the contact interactions would change observables even in 
the chiral limit. However, we have checked that the corresponding effects in observables are small and 
will not discuss this issue in what follows. The uncertainty in the LEC d\Q turns out to be much more 
important for our analysis. Contrary to the previously discussed case with the LEC d±s, changing 
d\Q result in changes of the strength of the OPE in the chiral limit (and, of course, for other values 
of M w / M w ) but does not affect it at the physical point where M n = M w . In what follows we will 
demonstrate how the uncertainty in diQ shows up in various observables. 

3.2 Phase shifts 

After the remarks of the preceding section we are now in the position to present results. In figs. 7 
to 10 we show our predictions for various phase shifts in the chiral limit compared to their values in 
the real world. In order to estimate the uncertainty due to missing higher order terms we proceed 
in a usual way and vary the cut-off between 500 and 600 MeV, as it has been done in [4]. The only 
exception to this are the 1 So and 3 S± — S D± channels in the chiral limit, where the dark shaded bands 
in figs. 7,8, 9 are due to the uncertainty in the values of As,T- Varying the cut-off in that case is not 
necessary, since the corresponding uncertainty is much smaller than the one due to unknown Ds,t- 
For these channels we will take the cut-off A = 560 MeV. 

Both S-wave phase shifts in fig. 7 show a qualitatively similar behaviour to what is observed in the real 
world. In the 1 Sq channel the interaction in the chiral limit gets weaker and the phase shift reaches in 
the maximum about 30 — 50% of the observed value. The situation in the 3 Si channel is just opposite 
but the effect is smaller in magnitude. Note that the band for the 1 Sq phase shift is much wider, 
partially due to the larger uncertainty in Di § compared to the one in D3 $ 1 ■ 

Before considering higher partial waves some comments are in order. First of all we note that the 
form of the OPE in the chiral limit 

{V )cl oc K 6 - 1 ) 

* 19 In our earliar work [3] we used a value for the pion-nucleon coupling g^N, which is slightly smaller than the commonly 
accepted one g n N = 13.1 . . . 13.4. We have refitted the contact terms Cs,t, C\.....i using the value for g n N according to 
eq. (3.1). This leads to minor changes in the values of the LECs related to contact terms and does not visibly change 
observables. 
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Figure 7: S-wave phase shifts as functions of the momentum in the cm. system. 
Dark shaded bands correspond to our NLO predictions in the chiral limit, light 
shaded bands are NLO results for physical value of M n , filled circles are Nijmegen 
phase shifts [30]. 



leads to significant scattering at vanishing momenta for all partial waves, as pointed out in by Bulgac 
et al. [10]. It is amusing to note that although the interaction between pions and nucleons vanishes 
in the chiral limit at vanishing momenta, the interaction between nucleons via the exchange of pions, 
becomes, on the contrary, stronger as a consequence of the increased range of the interaction for 
M w — > 0.* 20 As a consequence of the vanishing pion mass the effective range expansion 



where k is the cm. momentum and I the angular momentum, does not exist any more. The partial 
wave amplitude A\{E), E = k 2 /m,N, describing elastic NN scattering has a left-hand cut starting at 
the branch point E = — M 2 / {Am^), which gives a maximal radius of convergence of the effective range 
expansion. Thus the domain of validity of the effective range expansion ceases to exist for M n — > 0. 
Projecting the OPE in the chiral limit, eq. (3.1), onto the states with definite I, s and j we note that 
all spin-singlet matrix elements vanish. Further, on-the-energy shell matrix elements in all other 
channels do not depend on the momenta. As a consequence, one expects that the phase shifts at low 
energy in the high partial waves, in which the amplitude is essentially given by the Born term and 
strongly dominated by the OPEP, behave like 6(k) oc k. 

Let us now comment on the P-waves which are, as already pointed out in the introduction, of a 
particular interest, because of the possible existence of bound states due to the strong OPE in the 
chiral limit. First, we stress that we are able to perform an accurate extrapolation in the pion mass 
for P- and higher partial waves since the complete M^-dependence is known. Our findings confirm 
the conclusions of [10] that no bound states exist in the P-waves. The phase shift in the 3 Po channel, 
which is most sensitive to the OPE and thus is an ideal candidate for the appearance of a bound state, 
is strongly enhanced compared to the physically relevant case and reaches a maximum of about 32°. 
The l P\ phase shift becomes small and is completely given in the chiral limit by the TPEP and the 
NLO contact term. Remarkably, one observes strong changes in the mixing angle e±, which is known 

#20 The OPEP behaves at large distances in the chiral limit as 1/r 3 . 



k 2l+1 cot 5i(k) = -- + n lL + vfk 4 + ... 
a, 2 



(3.2) 
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Figure 8: P-wave phase shifts. For notations see fig. 7. 



to be an observable that is very sensitive to small changes in certain parameters. It is comforting 
to note that the bands do not get significantly wider compared to the case with the physical value 
of the pion mass, M n = M n . This is a clear indication of the consistency and good quality of the 
extrapolation to M n = 0, since the cut-off independence of observables in the chiral limit is achieved 
by "running" of the LECs Cj(A) fixed at the physical value M T = M n . 

The situation with the D- and F- waves turns out to be similar to the case of the P-waves. The 
phases in most channels (apart from the spin-singlet ones) already follow their asymptotical behaviour 
discussed above, i.e. grow linearly with momentum k. Our last comment in this section is that the 
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Figure 9: D-wave phase shifts. For notations see fig. 7. 



partial wave decomposition becomes not a reasonable (from the practical point of view) tool for 
calculating observables (full scattering amplitude) if M n — » 0, since the expansion in the partial waves 
converges slowly in the presence of the long-range interaction oc 1/r 3 . 

3.3 Deuteron binding energy 

We now move on to discuss the deuteron binding energy Bj), which we consider to be the most 
interesting observable with respect to its M^-dependence, see fig. 11. As already pointed out in 
the introduction, the situation concerning the chiral limit behaviour of £>d according to the previous 
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Figure 10: F-wave phase shifts. For notations see fig. 7. 



analyses of Bulgac et al. [10] and Beane et al. [11, 14] can not be considered as settled. In fact, the 
value of Bj) is rather sensitive to the assumptions and approximations made as well as to certain 
parameters, like the pion-nucleon coupling constant. According to our complete NLO analysis, we 
arrive at a unique result: the deuteron is stronger bound in the chiral limit with the binding energy 

#g L = 9.6±1.9 + j^ MeV. (3.1) 

Here, the first error refers to the uncertainty in the value of D:j Sl , where the LEC d\Q is set to the 
average value d\Q = —1.23 GeV -2 from [39]. The second error shows the additional uncertainty when 
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diQ is varied in the range given in [39]. We decided not to add these uncertainties but rather give them 
separately since they are not completely uncorrelated. The value of is in remarkable agreement 
with the natural value one would expect to arise in QCD of about F%/m ~ 10 MeV, if one (naively) 
assumes for the binding energy to be of the order of the kinetic energy with the relevant momentum 
scale being p ~ F^.* 21 




Figure 11: Deuteron binding energy as a function of the pion mass. The shaded 
areas correspond to allowed values. The light shaded band corresponds to our main 
result with die = —1-23 GeV -2 and the uncertainty due to the unknown LECs As,T- 
The dark shaded band gives the uncertainty if, in addition to variation of D$Ti the 
LEC dig is varied in the range from dig = —0.91 GeV -2 to die = —1-76 GeV -2 
given in [39] . The heavy dot shows the binding energy for the physical value of the 
pion mass. 



We also performed extrapolation for larger values of M n , which might be of interest for later calcula- 
tions within lattice QCD. According to our analysis the deuteron does not exist anymore for M n > 200 
MeV. Note that the behaviour of Bd is qualitatively similar to what has been found in [11] by taking 
into account only the explicit M^-dependence. 

3.4 S-wave scattering lengths and matching with lattice QCD 

Our last topic is related to the S-wave scattering lengths a So and a iSl and the possibility of matching 
with the calculations using lattice QCD. In fig. 12 we show our predictions for the scattering lengths. 
First of all, it is interesting to see that the uncertainty due to the unknown constant D3 Sl does not 

* 21 More precisely, one should use in this formula F and m, but for an approximate estimate the differences to the 
physical values can be ignored. 
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Figure 12: The S-wave scattering lengths as functions of M n . The shaded areas 
represent the allowed values according to our analysis. The light shaded band corre- 
sponds to our main result with die = —1-23 GeV -2 and the uncertainty due to the 
unknown LECs Ds,t- The dark shaded band gives the uncertainty if, in addition to 
variation of Ds : t, the LEC die is varied in the range from die = —0.91 GeV~ 2 to 
die = —1-76 GeV -2 given in [39]. The heavy dots corresponds to the values in the 
real world. The triangles refer to lattice QCD results from [13]. 



show up strongly for a 1 for M n < M n , where we are able to obtain an accurate extrapolation. In 

3 Q 

the chiral limit the scattering length a™ 1 is still positive and smaller in magnitude compared to the 
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physically relevant case, which is in agreement with the larger value of the deuteron binding energy. 
Specifically, we get = 1.5 ± 0.4 fm, which is again of natural size. Here the first uncertainty 
refers to variation in D$ t T while the second one to the additional variation in diQ, as described above. 
The scattering length changes its sign at M n ~ 200 MeV and relaxes to natural values for larger M n . 
Note also that the relative contribution of the OPE decreases with increasing pion mass and thus the 
uncertainty due to the variation in d\Q becomes very small. On the contrary, the uncertainty in Dss 1 
is, of course, magnified as M w grows. 

In case of the a s ° the uncertainty in our extrapolation in the pion mass is visibly larger. This is 
consistent with the known fact that the OPE plays a less significant role in this channel, where the 
interaction is dominated by the shorter range terms. Another reason is that an extreme fine tuning 
leading to the very large value for a So takes place in the physically relevant case with = M w , where 
all parameters are fixed. The larger effect due to the uncertainty in the DiSn thus has to be expected. 
Nevertheless, our results concerning this channel are also quite interesting. In particular, it turns out 
that no bound state appears for pion mass smaller than M^, although in the real world for M w = M w 
the virtual state in this channel is almost bound. Around the chiral limit a So shows a qualitatively 
similar behaviour to a Sl and gets smaller in magnitude: = —4.1 ± 1.6 fm, which is still 

somewhat large compared to what would be expected from dimensional reasons (1 — 2fm). We are 
not able to obtain an accurate prediction for a So for M w > M w . In that case both scenarios with or 
without a bound state in this channel are possible. This uncertainty leads to a huge uncertainty with 
respect to the scattering length, which becomes infinite in the presence of a zero-energy bound state. 
This situation is depicted in the upper panel of fig. 12. The different signs of the scattering length are 
due to the fact that the two-nucleon system might be bound or unbound in this channel. It is more 
appropriate in such a case to look at the inverse scattering length, which simply crosses zero when 
one encounters a zero energy bound state. We show the inverse scattering lengths 1/a s ° and 1/a' Sl 
in fig. 13. Unfortunately, it is not yet possible to compare our extrapolation of the scattering lengths 
with the results available from the lattice calculation. The latter ones are currently restricted to values 
of M n larger than 550 MeV, which is far beyond the domain of applicability of the chiral expansion. 
Actually, our extrapolation is not trustable anymore for M n > 2M n . For example, the relative shift 
A of the ratio gA/F-w defined in eq. (2.86) reaches ~ 50% for M n ~ 250 MeV. Furthermore, the lattice 
calculation of ref. [13] is carried out in the quenched approximation, so that one, in principle, should 
perform extrapolation in the pion mass in the quenched chiral perturbation theory, see [44, 45] for 
more discussion. 

3.5 Comparison to earlier work 

It is interesting to compare our findings with the ones of refs. [11, 14], which are also based on 
effective field theory. There it was pointed out that the (natural) values of the unknown LECs Ds,t 
exist, which make the deuteron both bound or unbound in the chiral limit and which strongly change 
the S-wave scattering lengths. Thus no predictions for all these quantities could be made [14]. The 
basic differences in our analyses can be summarized as follows: 

1. While we take the chiral limit value for the strength of the OPE, gA/F v , from the CHPT analysis 
of the 7T7T, irN and ttttN systems, the authors of ref. [14] make use of the assumption that terms 
of the type mM^/A, where M w <C A (A ~ M p ) dominate the chiral expansion. This leads to 
a decrease of qa/F^ in the chiral limit, while the CHPT based analysis gives an increased value. 
This difference seems to be the most important one. 

2. While we take into account the complete TPEP, only its chiral limit value is considered in [14]. 
In addition, we also include renormalization of the short-range interactions, which leads to the 



31 



100 200 300 400 500 600 700 

M n [MeV] 



Figure 13: The S-wave inverse scattering lengths as functions of M n . The shaded 
areas represent the allowed values according to our analysis. The light shaded band 
corresponds to our main result with d\Q = —1.23 GeV™ 2 and the uncertainty due to 
the unknown LECs Ds,t- The dark shaded band gives the uncertainty if, in addition 
to variation of D$ t T, the LEC d\Q is varied in the range from die = —0.91 GeV -2 to 
d\Q = —1.76 GeV -2 given in [39]. The heavy dots corresponds to the values in the 
real world. The triangles refer to lattice QCD results from [13]. 



nontrivial logarithmic dependence on the pion mass. 
3. We regularize the Lippmann-Schwinger in the momentum space by a finite cut-off and do not 
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require the low-energy phase shifts to be completely but only approximately cut-off independent. 
The remaining cut-off dependence can be eliminated by inclusion of higher order terms. In [14] 
the Schrodinger equation in the coordinate space (with the finite short-distance cut-off) is 
solved. The cut-off independence of the observables is achieved at the expense of introducing 
short-range interactions in the D-wave, which are not present in the effective Lagrangian. 

Finally, we remark that as in ref.[10] we find no new bound states in the P-waves. 

4 Summary and conclusions 

In this paper we have calculated properties of the two-nucleon system in the chiral limit based on a 
chiral effective field theory. The results of our investigation can be summarized as follows: 

• Based upon the modified Weinberg power counting (as explained in [2]) we have shown how to 
perform a complete renormalization within the method of unitary transformation, including the 
mass and wave-function renormalization and how to deal with tadpole graphs. We found the 
renormalized expression for the OPEP within the projection formalism and demonstrated that 
it agrees precisely with the off-the-energy shell extension of the OPE amplitude obtained in the 
S-matrix method. 

• Based on the NLO potential, which includes the renormalized OPE, TPE contributions and 
contact terms, we have performed extrapolations in the pion mass (or, equivalently, quark mass) 
away from its physical value * 22 . The corresponding LECs have been taken from an investigation 
of 7T7T [6], ttN [25, 39] and NN systems [3]. The only uncertainties arise from the unknown 
couplings -Ds,T related to contact terms with one insertion ~ M 2 . and to the variation in the 
pion-nucleon LEC dig related to the Goldberger-Treiman discrepancy. For the dimensionless 
coefficients that parameterize the LECs -Ds,t and which are expected to be of order one, we 
considered values in the range from —3 to 3. The variation in d±Q is obtained from considering 
various irN phase shift analyses as input to analyze the process irN — > ttttN. 

• We have calculated the NN phase shifts in the chiral limit. In the 1 S'o and 3 S*i — 3 D\ channels 
we obtain predictions with an error governed by the uncertainty in the values of I?s,T- Both 
S-wave phase shifts look qualitatively similar to the physically realized case, although the 1 S'o 
phase shift is about 50 — 70% smaller in magnitude. The 3 Pq phase shift is enhanced but no 
bound state appears in that channel, as it is also the case for other P-and higher partial waves. 
Starting from D- and F-waves, the phase shifts nearly reach their asymptotic behaviour caused 
by the exchange of a massless pion and grow linearly with momentum, S(k) ~ k. 

• According to our analysis, the deuteron is significantly stronger bound in the chiral limit. The 
binding energy is close to the natural value of ~ 10 MeV one expects to arise in QCD. There is 
no bound state in the 1 Sq channel in the chiral limit. 

• In the chiral limit, the S-wave scattering lengths take smaller (in magnitude) and more natural 
values as compared to the real world. Furthermore, for pion masses significantly larger than the 
physical one we found negative values of natural size for the 3 Si scattering length, which seem 
to be consistent with the lattice calculation [13]. In the 1 S'o channel the uncertainty due to the 
unknown LEC Di s is much larger than in the triplet channel and we were not able to predict 
o^So f or j ar g e values of the pion mass. The scattering length takes both positive and negative 

#22 It certainly would be interesting to extend the calculation to NNLO. 
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values with the magnitude changing from being natural to infinitely large. We also consider 
inverse scattering lengths, which are more suitable for chiral extrapolation and for comparison 
with the lattice calculation. 

• We have shown that it is possible that Wigner symmetry is also exact in the chiral limit and 
that its breaking observed in nature is entirely due to the quark mass terms. This would explain 
why the LECs related to Wigner symmetry breaking turn out to be so much smaller than all 
other LECs parameterizing the short-distance part of the NN interaction. 

To conclude, although we did not find dramatic changes in the properties of the NN systems in the 
chiral limit, the physically realized value of M n turns out to be a quite specific one in the sense 
that it leads to unnaturally large scattering lengths in both S-waves and consequently to the small 
deuteron binding energy. In fact, nature seems to be more simple in a world with massless quarks 
since nuclear binding would be of its expected size and the S-wave scattering lengths would be natural. 
Thus, no fine tuning would be needed and the nuclear force problem would be amenable to a much 
simpler treatment than required in the real world. Of course, our investigation can not answer the 
central question why nature chooses such a fine-tuned scenario for the nuclear forces, but at least it 
is comforting to know that she seems to be more kind in the chiral limit. 
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